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and that ^-^0 wheri to -^ along either of two radii vectores interior to A. Then 
<f)~^0 uniformly ivhen ia-^0 in any manner tmthin any region is! interior, to A. 
6.3. Similar arguments, based upon Lemma 5, lead to the conclusion that 
in MonteVs theorem of 6.1, the hypothesis that ^ (to) is hounded onay be replaced 
hy the hypothesis that the integral (1.15) is convergent. A direct proof of this 
theorem is not without interest, but the actual result is included in the more 
general theorems referred to in 6.1.* 



On the Free Transverse Vibrations of a Uniform Circular Disc 
Clamped at its Centre ; and on the Effects of Rotation. 

By E. V. Southwell, M. A. 
(Communicated by Prof. H. Lamb, F.R.S. Received September 23, 1921.) 

Inteoduction and Summaky. 

1 1. A thin circular disc, entirely unconstrained, can vibrate transversely in 
modes which may be characterized by any number of nodal concentric circles, 
crossed by uniformly spaced nodal diameters. The natural frequencies, for a 
disc of uniform thickness, have been calculated by Kirchhoff',f and a recent 
investigation by Prof. Lamb and the present author^ has shown the amount 
by which they are increased when the disc rotates at any given speed about 
its axis. 

This extension of Kirchhoff's analysis was suggested by published accounts 
of failures in turbine discs, the circumstances of which indicated that they 
had their origin in forced vibrations, intensified by resonance. | It was 
realised that the results would be only indirectly applicable to design, mainly 
because in practice discs are constructed with curved profiles (the thickness 
increasing from the rim to the boss), but also because the calculated 

* For (ji (w) covers a finite area when to covers A, and there are therefore infinitely 
many values which (j) (w) does not take. 

t "Ueber die Schwingungen einer. elastischen Scheibe," 'CreJle's Journal,' vol. 40 
(1850), and ' Pogg. Ann.,' voL 81 (1850). The essentials of the analysis are i^eproduced in 
Lord Rayleigh's ' Theory of Sound,' vol. 1, §§ 218, 219. 

J * Boy. Soc. Proc.j' A, vol. 99, pp. 272-280 (1921). Throughout the present investiga- 
tion this paper will be quoted by the abbreviated title "Spinning Discs," and its 
equations by their number with an A affixed : thus equation (34a) below is equation (34) 
of the '* Spinning Discs '' paper. 

§ Cf, K. Baumann, " Some Recent Developments in Large Steam Turbine Practice," 
Journ. Inst. Elect. Eng.,' vol. 59, pp. 619, 620 (1921). 
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frequencies, in some instances, will be seriously affected by the " clamping " 
which is introduced when the disc is attached at its centre to a shaft and 
rotor of considerable inertia. 

Varying thickness, and the allied complication of loads imposed at the 
rim by the inertia of the blading, introduce difficulties of such magnitude 
that it will probably not be possible to calculate the effects of these factors 
with any great accuracy. But the effect of central clampiug may be studied 
in a general fashion, if we confine our attention, as before, to a disc of uniform 
thickness, and calculate, by a further extension of Kirchhoff's analysis, the 
increase in frequency which results when a constraint of this nature is intro- 
duced. The present paper contains the results of an investigation on these 
lines. Its scheme is generally similar to that of the joint paper to which 
reference has been made ; that is to say, exact calculations are given for the 
two limiting cases obtained by neglecting in turn the centrifugal and the 
flexaral stress-systems, and these are followed by an approximate treatment 
of the problem in its general form. 

The assumption made in regard to " central clamping " requires some 
explanation. It is evident that in practice any alteration in the dimensions 
of the boss will modify not only the boundary conditions imposed upon the 
transverse displacement, but also the centrifugal stress-system, which in 
part determines the differential equation governing this displacement at all 
points in the disc. Both of these effects will infiuence the natural frequencies 
of vibration, the first in a way which we may regard as direct, and the second 
indirectly. N^ow the centrifugal stress-system, in practice, will depend upon 
many factors (e.g., the diameter of the boss, the diameter of the shaft, 
the stresses imposed by the process of ''shrinking on," the weight of the 
blading, etc.), and it is not possible to make any one assumption in regard to 
its distribution which will satisfactorily represent all cases. For the purpose 
of the present paper, which is intended to throw light on the general 
principles involved, it has therefore been thought best to concentrate atten- 
tion upon the direct consequences of clamping, and to assume that the clamps 
prevent transverse displacement without altering the centrifugal system, so 
that this system is such as obtains in a disc which is continuous up to the 
centre. 

A constraint of this nature is easily visualized {e.g., clamps which are 
brought into contact with the rotating disc after it has taken up the strains 
appropriate to its speed of rotation), and it is easy to reproduce in experiments. 
The approximate methods given in the concluding section of this paper 
enable us to deal with any other system of centrifugal stresses which may be 
specified, and in our exact calculations (§ 3) for the spinning membrane (or 
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disc of negligible flexural rigidity), by retaining unchanged the differential 
equation used previously, we are led to the interesting conclusion that in this 
limiting case clamping has no direct effect upon the frequencies, although it 
alters the modes by introducing terms which would become infinite at the 
centre ; it follows that the only way in which clamps can modify the 
frequencies is indirectly, by changing the centrifugal stress-system. 

The other limiting case, of the non-rotating disc, is discussed in § 2. 
Figs. 3-6 show how the frequency varies with the radius of clamping in the 
gravest modes having 0, 1, 2 or 3 nodal diameters, and some of the higher 
frequencies have been calculated for the special ease in which the radius of the 
clamping circle is infinitesimal. It was to be expected that a constraint of 
this pai^icular type would involve no alteration in frequency for modes having 
two or more nodal diameters, since in such modes both the displacement and 
the slope vanish at the centre even when the disc is completely free. It is 
also found, somewhat unexpectedly, to be without effect in modes having one 
nodal diameter,* although the frequency in the gravest mode rises very 
rapidly with an increase in the radius of the clamping circle ; but for the 
" symmetrical " modes, in which there are no nodal diameters, an entirely new 
series of frequencies is obtained when the centre of the disc is clamped. 

In the general case, when the flexural and centrifugal systems are both 
operative, the conclusion is reached that clamping confined to a small circle 
has a negligible influence upon the gravest frequencies in modes which have 
two or more nodal diameters, and that sufficiently exact estimates of these 
frequencies can be obtained, as explained in the '' Spinning Discs " paper, 
from the relation 

^2 :-: 2J{^ +P2^, (34a) 

where pi and p2 are the gravest values of p found by neglecting (1) the flexural 
rigidity and (2) the rotation. This formula gives an estimate which is of the 
nature of a loiver limit to the true value. It also yields reasonably accurate 
results as applied to modes having one nodal diameter, if p2 is given the 
value appropriate to the size of clamping circle considered. But in the 
symmetrical modes the gravest frequency found by neglecting the flexural 
rigidity is zero, the effect of the central constraint being nil] hence, the 
formula (34a), as applied to this case, would make no allowance for the effect 
of rotation. It is evident that the flexural rigidity, in a disc of which the 
centre is clamped, must in practice prevent the occurrence of modes in which 
the displacement, or its slope, is discontinuous at the centre, and will in this 
w^ay bring the centrifugal system into operation. Special calculations are 
therefore required for the symmetrical modes, and the gravest frequency has 

^ In these, dw/cr has finite vahies at the centre when the disc is completely free. 
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been investigated (Section lib, §5) by means of Eayleigh's method of an 
assumed type,* in which the potential and kinetic energies are calculated on 
the basis of an assumed form of w, and the frequencies determined from the 
condition that their sum must be constant in the vibration. It is known that 
this method gives an upper limit to the true value, and hence it may 
also be applied to modes having nodal diameters, as a convenient means 
of estimating the limits of possible error involved in the use of the 
formula (34a). 

Although it still has no direct bearing on the occurrences in actual turbine 
discs, it is evident that our analysis has been made more representative of 
practical conditions by the assumption of central clamping. The new- 
frequencies found for the symmetrical modes are of special interest, since 
it appears that vibrations of this type have been a source of trouble in 
practice.f 

Exception may perhaps be taken to the fact that in this, as in the former 
paper, certain factors which would require consideration in a strict theory of 
flexure have been ignored: no account has been taken of the additional 
deflections which accompany the transverse shear stresses, or of the terms 
which depend on " rotatory inertia." Both of these factors increase in 
importance with the thickness of the disc, and occur whether the disc rotates 
or not. In the latter case they are generally neglected, and are in fact quite 
unimportant except in the higher modes of vibration, where the disc is 
divided by nodal lines into a large number of segments. When the disc 
rotates, the terms which depend upon- '' rotatory inertia " are modified by the 
addition of terms of "gyroscopic" type, containing the angular velocity to as 
a factor. It is not difficult to show that these cancel out of the general 
differential equation of transverse vibration, but they appear in the boundary 
conditions. On the other hand, in the problems of this and the former paper 
they prove to be of the same order of magnitude as the other '' rotatory 
inertia " terms (of ordinary type), and may therefore be equally neglected. 

The author desires to acknowledge his indebtedness to Miss B. S. Gough and 
Miss S. W. Skan, of the Aeronautics Department, National Physical Laboratory, 
for the valuable assistance they have given him in the detailed calculations of 
this paper ; also to Miss C. G. C. Phillips, who prepared the diagrams for repro- 
duction. 



■^ 'Theory of Sound,' vol. 1, §§ 88, 89, 
+ (Jf. Baumann, loc. cit. 
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I. Feei Transyersi Vibrations of a Non-Rotating Circular Disc 

Clamped at its Centre. 

§ 2. Following Kirchlioff and Eayleigh, we shall neglect the effects of 
" rotatory inertia " and of the additional deflection caused by shear. Then the 
differential equation which governs the transverse displacement (w) of the 
middle-surface may be written, in polar co-ordinates, as follows^: — 



d'w 






E/i2 



di;' Z{l-a') 



3,.2 r 9r 1^ o&^A 



0. 



(1) 



In this equation, p denotes the density of the material, E Young's 
Modulus, and a- Poisson's Eatio, while li is the half -thickness of the plate. 
The boundary conditions of the present problem are 
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at the free edge (r = a), and 

3if; ^ )> at the circumference of the clamping circle (r = &). 



or 



0, 



(3) 



Assuming a solution of the type 

%v = W sin s {d— 6q) sinp {t—U), 
and substituting in (1), we obtain the equation 



(4) 



fP 1 d ^s 



2 



W-~/c^W = 0, 



where 



_dT^ r dr r^^ 



fC^ 



Wv" 



(5) 



(6) 



The complete expression for W thus consists of the two solutions of the 

equation 

r ^2 1 ^7 C.2 n 

W = 0, (7) 



'dr^ r dr r^ 



together with the two solutions of the equation 



dr'^ r dr 



s' 



.2 



fC 



2 



W = 0. 



(8) 



■^ Of. Bayleigh, * Theory of Sound,' voL 1, §218, or Love, * Theory of Elasticity,' 
§ 314(d), 
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The solutions of (7) are the Bessel functions denoted by Js(/cr) and Ys{k^t), 
and the solutions of (8) are the modified Bessel functions denoted by ls{fcr) 
and Ks {tcr). Thus we have 

W = a J. {kt) + ^Y, {r.r) + 7I. i^r) + SK, {kt\ (9) 

where a, /3, 7 and S are arbitrary coefficients. 

Then if we employ the recurrence formulae for Bessel functions, viz., 



n 



•^ 



^n {^) = - JnC^)-Jn+l(^), 






n 



y 



(10) 



In'(^) = - I«(^)4-I« + i(^), , 
Z J 



and the similar formulae for Yn (^) and K^' (^) respectively,^^ the edge 
conditions (2) may be written in the formf 



K? % Ot 



^ 



J, (««) — (1— cr) \ '!-£-— -2 Jj (/<;«) + -_J,+i(A:a) 
Y. («a) _(l_a) \'-^~^ Y,(«a) + i- Y,+i («a) 



L fc^a^ 



-7 



--S 
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* WMttaker and Watson, * Modern Analysis,' §§ 17.21, 17.61, 17.7, and 17.71. 

t The coefficients k and jc^ are important, since the equations may be satisfied solely by 
reason of k being zero: a being necessarily finite, equation (12) has been multiplied 
throughout by this quantity. 
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whilst the conditions (3) at the circumference of the clamping circle become 



and 
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(13) 



0. (14) 



Eliminating the coefficients a, /8, % S from (11) — (14), we obtain a deter- 
minantal equation which, for any required value of the ratio b/a, may be 
solved by trial,* assumed values of fca being adopted, and the resulting values 
of the determinant plotted against these. Table I gives results obtained in 
this way in respect of the gravest frequencies for modes in which s has the 
values 0, 1, 2 and 3 : cr has been taken to be 0*3. 



Table 1. 



s ~ 


= 0. 
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2-50 
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8-00 
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0-771 


8-00 


0-81 


10-00 










0-827 


11 -00 











In the limiting case, when the radius of the clamping circle is indefinitely 
reduced, the analysis may be considerably simplified. For if s has a value 
greater than 1, the coefficients of a and ^ in (13) and (14) vanish when 
& = 0, and these relations will accordingly be satisfied by making 7 and S 
zero, without restriction of a or /3. The frequencies are thus found by 
eliminating a and ^ between (11) and (12), which now reduce to the 
corresponding equations of Kirchhoff for a completely free disc. 



* Tables of the functions J^{z) and ln{z) (n integral) are given by Gray and Matthews, 
* Treatise on Bessel Fimctions,' pp. 247-288. The Neumann functions Yq (z) and Yj (z), 
as here defined, have been tabulated by the British Association, ' Eeport for 1914,' 
pp. 75-82, and the functions Kq (s) and K^ (z) by W. S. Aldis, ' Eoy. Soc. Proc.,' A, vol. 64, 
pp. 203-223 (1899). 
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When s ha,s the values 1 or 0, we may substitute in (13) and (14) from the 



ascending series for the Bessel functions. We have 



Jo(^') — 1 — -J- + ..., 



J, M = ~ _- + 



(V 



z 



3 



U(z) = 1 + J + ..., h(^) == '^'^iq"^'"' 



z^ 



Yo(z) == Jo{z)logz + -^ + ..., 






Yi(,^) = Ji{z)logz~~-^~-j + z^,, 

Z 4fc 



y 



(15) 



Ko(^) 



z' 



lo(^) {log ;<J — E}4- - + ..., 



Ki(2;) = — Ii(:^){log2;— E} 



1 . z 



-f- ■- 4" • • • 



;s 



where 



iSi. 



E = log 2 — 7 (7 being Euler's constant), 



y 



Hence (observing that the limit of ;c;log2;, when z-^O, is zero) we see that 
the signijficant terms in (13) and (14), when s = 1, are 

1 



fcb 



(/3 + S)=:0, 



1 



and 



2 • ^\~2~^«252 + 4, 



r./'l0S«& — E 

H— 2 



1 1\ 
K'^b'' 4/ 



^ (16) 



V 



respectively. It follows that (/3 + S) must vanish, and that both ^ and S must 
be infinitesimal in comparison with a and 7, although their contribution to 
the value of the slope is finite at the centre. Hence, if /ca is finite, the con- 
tributions of y8 and S to (11) and (12) will be infinitesimal, and the finite 
frequencies of vibration in modes characterized by one nodal diameter can be 
found from a relation identical with that which obtains for a completely free 
disc ; that is to say, these frequencies also are unaffected by clamping wdiich 
is entirely confined to the centre. 

The gravest frequency for the completely free disc, in modes characterized 
by one nodal diameter, is zero, the "vibration" consisting simply in a rigid- 
body rotation of the disc about a diameter. We should expect that this zero 
frequency will become finite when the slope of the disc is constrained at the 
centre; but it is easily verified that equations (11) — (14), when 5=1, are 
satisfied if we write 

/ca z=: fcb =^ 0, 

* This significance for E (the symbol employed in the tables referred to) continues 
down to equation (24), after which E is again used to denote Young's Modulus. 
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so that a zero frequency of vibration is still possible in the limiting case, when 

the radius of the clamping circle is indefinitely reduced. Alternatively, the 

existence of a mode having zero frequency can be established by verifying that 

the expression 

tv = {er -f ^r log r) sin (0 — ^o), 0-^) 

where e and f are arbitrary constants, is a solution of equation (1), which 

will also satisfy (2) and (3), provided that ^ is infinitesimal, but related to e 

by the condition 

€+f[logr + l],.o-0. (18) 
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Fig. 1. 



On the other hand, when the results (Table I) obtained for finite values of 
h/a are plotted^ as in the curve AB of fig. 1, kcc appears to be tending to a 
finite value as h/a-*0. It thus appears that the frequency must rise very 
rapidly as the radius of the clamping circle increases from a zero value, and 
the conclusion may be confirmed by expanding the determinantal equation 
which is obtained by eliminating a, 13, 7, S from (11)~(14), on the assumption 
that Kb is very small. When 5=1, the determinantal equation is easily 
thrown into the form* ' 

A^ AJi(A)-7;J2(A), AYi(A)-7;Y2(A), AIi(A)--7;l2(A), AKi (A)-7?K2(A), 



AJi (A) - 


-V^2(A), 


AYi(A)- 


-7?Ya (A), 


AIi(A)- 


-vh{A), 


AKi (A) 


J2 (A), 




Y2 (A), 




-12(A), 




-K2(A), 


BJi (B), 




BYi (B), 




-BTi(B), 




-BKi(B), 


Jo (B), 




Yo (B), 




-Io(B), 




-Ko(B), 



0, 



(19) 



■^ a being finite, it is legitimate to multiply the first two rows by this quantity. The 
complete equation contains an additional factor 1/6, which has been cancelled out of 
equation (19), since it cannot be zero. 
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in which A and B are written for Ka and icb respectively, and t] is substituted 
(for brevity) in place of the quantity (1-— o"). If we expand the determinant, 
substituting from (15) for the Bessel functions of argument B, the positive 
powers of B are found to be associated with powers of A not lower than 
the fourth, whilst the terms which are important when B ~> may be obtained 
from the deter minantal equation 

A^ AJi(A)-7;J2(A), AYi(A)~7;Y2(A), AIi (A) -7712(A), AKi(A)-7;K2(A), 

J2(A), Y2(A), "-12(A), "K2(A), 

0, -1, 0, 1, 

1, logB, -^1, (logB-E), 

= 0. (20) 

Hence we find that the complete equation (19) may be written in the form 
(positive powers of B) x (terms in A'^ and higher powers of A) 

— [(3 + a) A^ 4- (higher powers of A)] 

A« 



(21ogB-E) 



'^.2 



(3 + (t) -f (higher powers of A) 



= 0, 



.0. 
and when A -»► the significant terms in this equation are 

-(3 + <T)A4/l+^'logB)=:0. (21) 

The solution A = (for all values of B) corresponds to the configuration 
of equilibrium, with which we are not at present concerned ; the alternative 
solution may be written (since A log xi -> with A) in the form 

A4 = -r^v where B = hi a, (22) 

log E 

Hence, A -* when jR -> (the result found previously), but 

B V \BHogB)' 

so that — /or --— ) -* 00 as ^ -^ 0. 

B \ clBI 

The diagram shown in fig. 1 connects values of A (or tea) with values of B 
(or Ijci) ; AB is a curve drawn through points which have been plotted 
from Table I, and OC is a curve obtained from the approximate determinantal 
equation (20), tables being employed to find the exact values of the Bessel 
functions in the first and second rows. In this connection, it is convenient 
to employ the substitutions 

AJ2 (A) = 2 Ji (A)-- A Jo (A), 
AY2 (A) = 2Yi (A) - AYo (A), 
^Al2(A) = 2Ii(A)~-AIo(A), 
-AK2(A) = 2Ki(A)-~.AKo(A). 
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From the fact that the two curves very nearly overlap, it is evident that (20) 
approximates closely to the correct equation over the range considered ; a 
sufficiently exact transition from one curve to the other may be made by eye. 
This completes the discussion of clamping for modes in which s = 1. When 
s = 0, the significant terms in (13) and (14) are 

ot + 7 + (i8--S)log«:& + SE = " 
and fc^=0. ^ <^^^ 



Kb 
Both relations are satisfied if we write 



/3 = S = — -(a + 7X 



(24) 



and (11) and (12) can then be expressed in terms of u and 7 only. Eliminating 
a and 7, we obtaia an equation in kci which may be solved, as before, by trial. 
The results depend to some extent upon the value of a, which throughout 
this paper has been taken to be 0'3 (a representative value for steel) ; for this 
value, the first four calculated roots for Ka are 

Kct = 1-937, 4-573, 779, 10-94, 



and the corresponding frequencies can be obtained from (6) in the form 



F 



pa' 



X (5-156, 160-2, 1349, 5250), 



(25) 



where E is Youns^'s Modulus for the material considered. 



'D 



Fig. 2 illustrates the nature of the two gravest modes. 




Fia. 2. 

The variation of the frequency with the radius of the clamping circle, as 
found from the foregoing analysis, is exhibited in figs. 3-6, by means of 

curves which connect JPa/ (^^/(^ non-dimensional quantity proportional 

to the frequency) with the ratio B(=: h/a) ; these show, as we should expect 
from physical considerations, that the effects of clamping are unimportant for 
small values of B, unless s = 1, when Ithey increase very rapidly with the 
radius of the clamping circle. 
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II. The Effects of Eotation. 

(a) Fleamral Rigidity Negleetecl, 

I 3. As in the paper on " Spinning Discs/'* we begin by considering the 
extreme case, in which the disc is so thin, and the rotation so rapid, that the 

"^ Loc, cit,y § 2. On the assumption stated in the introdiictorj section of this- paper, 
the centrifugal stress-system, and lience the equation of transverse motion, will be 
unaffected by clamping. 
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flexural forces are negligible. Under these conditions, the equation of trans- 
verse motion may be written in the form* 

A(a^-O(|H+l^)_2Ar|-"+f4+«^B-£?^>. = 0, (5a) 

\dr^ r ar / or \6>^ r^ / 

where A = i (3 + o-), B = -| (1 + Sa). (2a) 

8 8 



If we write 



£--f s^B = m(m + 2)A, (8a) 



and (^ ■--- — 5j 

equation (5a) may be reduced to the form 

CW~\ . f / . cx\ s 

Again, if we write 



9 r P.«"^ r o2i 



4±L J^(l«.^)^if L + J rii(m + 2)--- }> w = 0. (26) 



in (26), we obtain, as the differential equation to be satisfied by y, 

4^(l-.^)^4-4{l + s-~(s + 2)4^+{m(m + 2)~5(5-|-2)}2/=: 0, 
and this may be written in the form 



CiX vO V J- *""* X ) iX/JO t^ \^ X """' X I 

where a = | (s— m), ^8 = H^ + ^+2), ry = s-f 1. (12a) 

Equation (27) is hypergeometric. The solution finite for a? = is 



1-7 1.2.7(7 + 1) 



y (28) 



^ 



where C is an arbitrary constant, 
or fi =CF(a, 0,%m) 

(in Gauss' notation). This is the solution which was discussed in the paper 

on "Spinning Discs," and the natural frequencies were obtained from the 

known theorem that, since 

ry = a + /3, (13a) 

the expansion of E (a, ^, 7, x) becomes logarithmically infinite at the edge 
(^ = 1)^ unless it tormina tes.f Hence we must have 

m = s + 2n, (14a) 

where ^ is a positive integer. 

•* The suffix previously attached to p is for the present omitted, 
t Forsyth, 'Differential Equations,' § 113, 

VOL. CI. — A. , L 
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To take account of constraints which prevent transverse displacement 
along some circle {x = constant), we must consider the second sohition of 
(27). When s is finite; this may be written in the form* 

where a — a + 1— 7 = — |(5 + m), 

^' = /3 + l^ry==i(m + 2--s), j 

ty' = 2 — ly = 1 — S, J 

whilst if s = (so that 7=1) the second solution isf 

2^ 



y (29) 



y, = C [f («, ^, 1, X) \o^x+ 1^ {~ + \- 



r" 



«(« + l)i3(/3 + l) /l 1 1 _J 2_2' 

. 1^.22 V« a + 1 /8 /3+1 1 2/ 



,(30) 



C denoting in each instance an arbitrary constant. Now, since 

7' = ^' + /3', 

the expansion of F (a', /S'/^', x), in (29), is divergent unless it terminates. But 
it is easily seen, from (28), that if F (a, /3, 7, x) terminates with the term in 
x"", so that (14) A is satisfied, then F {a\ /3', 7', a?) will terminate with the term 
in a?^+^so that the solution given by (29) consists of a finite number of terms, 
beginning with a term in x"^ and ending with a term in x^. 

Again, it is evident that the quantity in twisted brackets on the right of (30) 
will terminate under the same conditions as F (a, yS, 7, x), and at a term of the 
same order in x ; moreover, the factor F (a, j3, 1, x) in the first term is the 
solution found in (28). Hence, for all values of s, by suitably choosing the 
ratio of C to G in the complete solution 

we can make y (and hence w) vanish at any given radius ; and although the 
mode will he altered, the frequency of vibration will he unchanged. When the 
constraint is applied at the centre, so that the radius at which w vanishes is 
indefinitely reduced, it is clear that the ratio of C to will be infinitesimal, 
so that the mode is unaffected by the constraint, except at the centre : at the 
centre it will be affected only if 5 = 0, since for other values of s the central 
displacement is zero in the first solution (28). 

(b) Flexural and Centrifugal Stresses both Operative, 

§4« When the centrifugal tensions and the flexural stiffness are both 
operative, the difficulties of an exact solution are much increased, since the 

^ Forsyth, op. cit., § 115. 

t Forsyth, op, cit., p. 251, Ex. 1. 
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modes as well as the frequencies of vibration depend upon the relative 
importance of the two systems. In regard to the gravest frequency corres- 
ponding to any given value of s, it was shown in the paper on " Spinning 
Discs " ^ that a loiver limit to the frequency* p can be obtained from the 

equation 

^2 -^ ^^2 _|.p^2^ (34a) 

where _pi and f2 are the values of p obtained by neglecting the flexural 
rigidity, and the effects of rotation, respectively. 

The discussion just given shows that the values of pi^ are unaffected by 
clamping. Thus we may write, as in the earlier paper, 



Fi^ = Xm^, 



(20a) 



where X has values which may be found from the relations (2a), (8a) and 
(14a), given above. The following table (taken from the earlier paper, and 
amplified by the inclusion of a column for modes in which s = 0) gives 
numerical values for X on the hypothesis that a = 0*3, the value assumed 
previously :— 

Table II. — Values of X in the formula p-^ = Xo)^. 
(N.B. — These values are not affected by clamping.) 





^ = 0. 


s ^1. 


.==2. 


5 = 3. 


% ^ 





1 


2-35 


4-05 


n = 1 


3-3 


5-95 


8-95 


12-3 


n = 2 


9-9 


14-2 


18-85 


23 -85 


*i = 3 


19-8 


25-75 


32-05 


38-7 



The results of Section I show that the radius of clamping has an influence 
upon the value of p2^ but that when this radius is indefinitely reduced the 
only frequencies affected by clamping are those corresponding to the modes in 
which s = 0. Hence, writing 



V2' 



/^ 



pa^ 



(2lA) 



as in the " Spinning Discs " paper, we have the following table of values for 
/fc, of which only the second column (s = 0, centre clamped) is new : a is again 
assumed to have the value 0*3. 



"^ Log. eitf § 4, 
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Table III. — Values of jjl in the formula pi = jjl 






(N".B. — Clamping is confined to the centre of the Disc.) 





s = 


-- 0. 


s - 1. 
(Tliese T 


,9 = 2. 


6' - 3. 


Centre 
free. 


Centre 
clamped. 


aiues are not affected by 
sntral clamping.) 


% == 

n — 1 
n ^2 
n = 3 



29-67 

577 
2983 


6-16 
160-2 
1349 
5250 


0^^ 
164-3 
1222 
5266 


10-61 
443-9 
2624 

8574 


56-87 
1040 
4285 
13287 



* But cf. the discussion of this case given in §2. When the radius of clamping is finite, 
though very small, 1 '06 may he taken as a representative value of fx. 

By means of these tables (or of Table II and figs. 3 — 6, when the radius of 
clamping is finite) the lower limit given by (34a) can be quickly calculated. 
But it is evident that no allowance will be made in this way for the effect 
of rotation in regard to the modes in which s — : these, in a disc of which 
the centre is clamped, must be conditioned in the main by the flexural system, 
and will necessarily be such as will briug the centrifugal system into 
operation. 

§ 5. In the " Spinning Discs " paper,* closely approximate results were 
obtained by means of Kayleigh's method of an assumed type,t in which the 
potential and kinetic energies are calculated on the basis of an assumed 
expression for w, and the frequencies determined from the condition that their 
sum remains constant in the vibration. A frequency thus calculated is 
stationary for slight variations of the type, and a good approximation can, 
therefore, be obtained in this way if the type be suitably chosen. Moreover, in 
the case of the gravest mode of all, it is known that the frequency will be 
overestimated, and in the present problem this statement applies also to the 
gravest mode corresponding to any definite value of s. Prof. Lamb has made 
use of this fact in devising a method whereby the accuracy of Eayleigh's method 
may be improved ; it consists in the inclusion of a variable parameter in the 
assumed expression for w, the value of this parameter being afterwards 
adjusted so as to make the calculated frequency a minimum. J Using this 

■^ Loc, cit,^ p. 276. 

t * Theory of Sound,' voL 1, §§ 88, 89. 

t Cf. his paper " On the Vibrations of an Elastic Plate in Contact with Water," * Eoy. 
Soc. Proc.,' A, vol. 98, pp. 205-216 (1920), which appears to be the earliest illustration 
of the method. 
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device, the gravest frequency of vibration in a mode having two nodal 
diameters was estimated, for a non-rotating disc, with an error less than 
I per cent.* 

Applying these methods to the determination of the gravest frequency in 
modes for which s = 0, we have as expressions for the potential and kinetic 
energies, in polar oo-ordinatesf : — 

Vi = Potential Energy of the Centrifugal Forces, 

= 2-rrh^vi^^jrdr, (31) 

where P = A(a^~7^^) /ow^ (1a) 

and A is given by equation (2a) above : 
V2 = Potential Energy of Flexure, 






_dr^ r dr 



2 2(l — fT)dvj d^w 



r 



a;: • 3-2 f ^^^ ; (32) 



T = Kinetic Energy of the Motion, 

= 2,rph (^] TclT. (33) 



In the application of Eayleigh's method, an assumed form for w is substituted 
in these expressions, and we obtain an estimate of the frequency from the 

energy equation 

V14-V2 + T = constant. (34) 

We begin with the extreme case in which m is zero, since accurate results 
are here available for comparison. The simplest expression for w which 
satisfies the conditions that to and dioldr vanish at the centre is 



ni^ 



w = k ~ sin^2^. (35) 



a 



where k is an arbitrary constant. Substituting this expression in (32) and 
(33), we have 

y 

T = -. k^p2^ pha^G08^p2t, 



3 



and from (34) we obtain the value 



2 8 EA2 



1 — o" pa'^ 
whence fi = 11*429, when a = 0*3. (36) 

* Vide § 5 of the " Spinning Discs " paper, equation (52). 
t Cf, equations (25)~(29) of the " Spinning Discs " paper. 
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Comparing this result with the correct figure (5'16) given in Tahle III, we 
see that it is over 100 per cent, in error. To obtain a closer approximation, 
we may employ the artifice of Prof. Lamb, and assume, in place of (35), that 



Q.2 / T^\ 

■IV = « -J 1 4- /9 — , sin ■■p2t, 



(37) 



u and ^ being arbitrary constants. Then we have 

V2 = 



3^f^|j5{(l + ^)(l + 4^) + |(5 + 3^)/3^ h sin^M 



'^ 



3 2 5 



T = lT(^^i plhO? A ~^j\-^-\-~ \ C0S^j?2<^, 



)■ 



J 



and (34) gives the expression 



n 



^ 



8 EA^ f(l + a-)(l + 4/3) + 4/3 (5 + 3a-) Q^' 
3(1-0-2) /)66*L 



1/3 + /3/2 -1-/375 
Thus we have, in equation (21a) above, the expression 

L4-2M/3 + N"^2 

where, if o- ~ 0*3, 



(38) 





'" I + 2m/3 + n^^ ' 






T - 


3-81 


I - 


1 -] 
3' 


\j — 


— '■ — ' " — -...■.—. __ ^J (_jx, 

3(1-^) 


M = 


= ^ = '■''■ 


m 


.1 y 

' i' 


.N = 


._?2^^,,,„_ 


n — 


1 



(39) 



(40) 



or 



The stationary values of fi are given by the quadratic 

1 



0, 



240 



/^2--.4"64y^ + 29-79 = 0, 



whence the gravest value of j^2^ in (38), is given by 

p2^ = 6*457 -—^, approximately. 



(41 



a result which is still 25 per cent, in error. 

It is evident that the additional boundary conditions at the centre of the 
disc increase the difficulty of finding a suitable assumed form for tv. The 
value of ^ which corresponds to (41) is — 0-275 approximately, and when 
this is substituted in (37) the resulting form for to is found to„ differ materially 
from the correct form shown in fig. 2. After several failures to guess a better 
form., the plan was adopted of calculating the deflection of the disc when 
clamped at its centre and subjected to a loading varying in intensity as 
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some power of the radial distance, the edgej being assumed to he free from 
shearing forces, hut not necessarily from stress-couples. This deflection, 
multiplied hy sin ^2^, was taken as the assumed form of ^v, and the frequency 
was then calculated hy Eayleigh's method : finally, the index of r in the 
assumed intensity of loading was chosen so as to make the calculated 
frequency a minimum. In this way, the advantages of Prof. Lamb's device 
are retained, whilst the assumed mode, independently of the value given to 
the variable parameter, satisfies three out of ^the four boundary conditions. 
It is easily verified that the expression 



W = C -<^ r? - 



— a^ ^r^ log - 
4: a 



(42) 



where C is an arbitrary constant, gives the displacement corresponding to a 
loading of intensity 

(43) 



D 



1 d_ ^, dy 

T dr dr. 



W = GI)q^{q-2yf 



_0\2^.g-4 



2 EA^ 
where D =: - ^ ; it may also be shown to satisfy the second of the edge 

O JL — or 

conditions (2), together with both of the clamping conditions (3), provided 
that q^l and q ^ 2. 

For a value of ^ close to 3 it will also satisfy the first of (2); but' it is 
preferable to leave q indeterminate for the present, and to choose it 
subsequently so as to make the calculated frequency a minimum. We 
find that 



a 



r 







-dm,idw~^^^ 

— — — — —J— — — I 

dr^ r dr 



dr = CY^2(g-l) 



1 __2(q-^l)(q^2)_,(q-^2y 



2-i 



L2(^-l) 



r 



-f 



4 



and 



■« dW d^W 



dr= GYa^<^-^^ 



dr ' di'^ 

whence the potential energy of flexure in the mode 

■w = Wsin^i^ 



1 g(g-2) ,g'^(g-2yi 



32 



(44) 



IS 



3(l-o-2) 



(y\fa? <*~^^ sn)?pt 



.2(^-1) 
-(1-a) 
The kinetic energy of the motion is 

fa 
ry^^dr, 


= 27rphp^ GV ^^ "^ ^) cos^ pt 



2(s-l){q-2) + ti2^' 



g(g-2) I r/(g-2)n 
2 16 J _ 



(45) 



1 



qHq-2) q-^{q-2y 
2{q + l)' 2{q + 4:y 16. 108 J 






(46) 
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The expression for the frequency is therefore given by equation (21a) 
above, if 



/^= o 



r 



3(1-0-2) 



X 



5' 



9J. 



(^-1) 



afe-Dfe-aj+JM-'-d-,) (i-£fc|)+rtJj=IZ} 



1 q^{q^2 ) q^q-'^y 

2(^ + 1) "^2(^ + 4)2"^ 16 . 108 



(47) 



20'0 






o-o 




4'0 



-fin 



H 



Fig 7 



We have seen that the conditions to be satisfied at the centre require that 
q shall be greater than 1 and not equal to 2. Subject to this restriction, we 
may choose q so as to make ix, in (47), a minimum ; the required value is 
best determined by plotting fju against ^, as shown in fig. 7, where the 
important part of the curve is also plotted on a much enlarged scale. 
Taking a as equal to 0*3, it is found that the expression ^ has a minimum 
value of 5*22 when q = 2*89. This is an error of less than 1*2 per cent, in 
the value of fju, or less than 0*6 per cent, on the frequency, which is a very 
satisfactory result : but the rapid rise of /a, as q departs from the optimum 
value, serves to explain why such serious overestimates were obtained from 
what might have been considered reasonable assumptions for the form of %v. 

We proceed to consider the effects of rotation. Substituting for w, in 
(31 ), from (42) and (44), we have 



Vi = 2iTpm^A}hG\^a^^-^^d-D?pt 



1 q^{q-2) ^qHq-2f 



+ 



2 



L2^ 2(^ + 2)2 

1 ^ 

^ ^2(^ + 4)^ 



8.16 



q{q + 2)(q^2) 5qHq^2n 

i^+V) 2(^ + 4)2 3. 18. 16 J 
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The kinetic energy is given by (46), as before. Hence, by (34), the 
•expression for the frequency will be of the form (20a), where 

1 



\ = Ac/ 



^ (^-2) (r/ + 2q-A) ^ 7q^q--2)^^ 



2q(q~\-l) (q-^2f(q + 4) 



2 



16 . 27 . 8 



+ ^lill2I+2^(^-2)^ 



(48) 



2(^ + 1) 2(^ + 4)2 16 . 108 

\ as determined from this expression, may be plotted against q, in the 
manner of fig. 7. Its minimum value is 0*416, and occurs when 5' =3*295: 
like that of ytt, it is sharply defined, and the optimum value for q will clearly 
depend in the general case upon the relative importance of the flexural 
and centrifugal terms. But if we regard it as determined entirely by the 
flexural system, the value of X will be only slightly, if at all, in excess 
of 0*5 ; adopting this value as sufficiently accurate for practical purposes, 
we may calculate the gravest frequency in a symmetrical mode from the 
formula 



where 



2J^ = Xco^-i-fi 
0<X<0-5, 






y 



(49) 



and 5-16 < z^- < 5*22. 

A closer approximation could, of course, be obtained, in any particular 
instance, by using the curves for X and /x. Alternatively, by assuming the 
mode to be unchanged by rotation (so that fi is given by Table III), and 
employing (31), (33) and (34) for the calculation of X, we could obtain upper 
limits to the value of this quantity, for use in the general formula (49) : 
Table II, as we have seen, gives lower limits. The calculations would need to 
be done graphically. 
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